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^ • [Invited talk at TNT2005 (Trends in Nanotechnology), Physica status solidi (a) 203, 1160 (2006)] 

We discuss symmetry breaking in two-dimensional quantum dots resulting from strong interelectron re- 
pulsion relative to the zero-point kinetic energy associated with the confining potential. Such symmetry 
breaking leads to the emergence of crystalline arrangements of electrons in the dot. The so-called Wigner 
molecules form already at field-free conditions. The appearance of rotating Wigner molecules in circu- 
lar dots under high magnetic field, and their relation to magic angular momenta and quantum-Hall-effect 
, . fractional fillings is also discussed. Recent calculations for two electrons in an elliptic quantum dot, using 

' ^ ' exact diagonalization and an approximate generalized-Heitler-London treatment, show that the electrons 

C ' can localize and form a molecular dimer for screened interelectron repulsion. The calculated singlet-triplet 

O I splitting (J) as a function of the magnetic field (B) agrees with cotunneling measurements; its behavior 

^, . reflects the effective dissociation of the dimer for large B. Knowledge of the dot shape and of J{B) allows 

determination of two measures of entanglement (concurrence and von Neumann entropy for indistinguish- 
able fermions), whose behavior correlates also with the dissociation of the dimer The theoretical value for 
►^ I the concurrence at B = agrees with the experimental estimates. 

^nJ Copyright line will be provided by the publisher 
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^ ! 1 Symmetry breaking in quantum dots 

-(— > 



Two-dimensional (2D) quantum dots (QDs), created at semiconductor interfaces through the use of litho- 
graphic and gate-voltage techniques, with refined control of their size, shape, and number of electrons, are 
often referred to as "artificial atoms" LL.2.^- These systems which, with the use of applied magnetic 
fields, are expected to have future applications as nanoscale logic gates and switching devices, have been 
^ ■ in recent years the subject of significant theoretical and experimental research efforts. As indicated above, 

Q , certain analogies have been made between these man-made systems and their natural counterparts, sug- 

O ■ gesting that the physics of electrons in the former is similar to that underlying the traditional description 

of natural atoms - pertaining particularly to electronic shells and the Aufbau principle in atoms (where 
ij ■ electrons are taken to be moving in a spherically averaged effective central mean-field potential). 

rS I The above-mentioned analogy has been theoretically challenged recently |4j O on the basis of calcu- 

j^ ■ lations that showed evidence for formation, under favorable conditions (that are readily achieved in the 

laboratory), of "electron molecules," which are alternatively called Wigner molecules (WMs) after the 
physicist who predicted formation of electron crystals in extended systems (SJ. These spin-and-space 
(sS) unrestricted Hartree-Fock (UHF) calculations (denoted in the following as sS-UHF or simply UHF) 
of electrons confined in 2D QDs by a parabolic external potential led to the discovery of spontaneous 
symmetry breaking in QDs, manifested in the appearance of distinct interelectronic spatial (crystalline) 
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Fig. 1 The calculated spectrum of a two-electron parabolic quantum dot, with Rw ~ 200. The quantum numbers are 
(A'^, M, n, m) with A'^ corresponding to the number of radial nodes in the center of mass (CM) wavefunction, and M is 
the CM azimuthal quantum number. The integers n and rn are the corresponding quantum numbers for the electrons' 
relative motion (RM) and the total energy is given by ENM,nm ~ Ej^f + E'^"'(n, \m\). The spectrum may be 
summarized by the "spectral rule" given in the figure, with C = 0.037, the phonon stretching vibration huis = 3.50, 
and the phonon for the bending vibration coincides with that of the CM motion, i.e., huib = huio ~ 2. All energies are 
in units of hu)o/2, where luq is the parabolic confinement frequency. 



correlations (even in the absence of magnetic fields). Such symmetry breaking may indeed be expected 
to occur based on the interplay between the interelectron repulsion, Q, and the zero-point kinetic energy, 
K. It is customary to take Q — c^/kIq and K = Tllliq, where ^o = {h/m*uJoy^'^ is the spatial extent 
of an electron in the lowest state of the parabolic confinement; m* is the electron effective mass, k is the 
dielectric constant, and ljo is the frequency that characterizes the parabolic (harmonic) confining poten- 
tial. Thus, defining the Wigner parameter as Rw = Q/K, one may expect symmetry breaking to occur 
when the interelectron repulsion dominates, i.e., for Rw > 1. Under such circumstances, an appropriate 
solution of the Schrodinger equation necessitates consideration of wave functions with symmetries that 
are lower than that of the circularly symmetric QD Hamiltonian. Such solutions may be found through 
the use of the sS-UHF method, where all restrictions on the symmetries of the wave functions are lifted. 
On the other hand, when Rw < 1, no symmetry breaking is expected, and the sS-UHF solution collapses 
onto that obtained via the restricted Hartree-Fock (RHF) method, and the aforementioned circularly sym- 
metric "artificial-atom" analogy maintains. From the above we note that the state of the system may be 
controlled and varied through the choice of materials (i.e., k) and/or the strength of the confinement (ujq), 
since Rw oc l/{Ky/Lu^). 



1 . 1 Two-electron circular dots 

To illustrate the formation of "electron molecules," we show first exact results obtained for a two-electron 
QD, through separation of the center-of-mass and inter-electron relative-distance degrees of freedom [7J. 
The spectrum obtained for Rw = 200 (Fig. 1), exhibits features that are characteristic of a collective 
rovibrational dynamics, akin to that of a natural "rigid" triatomic molecule with an infinitely heavy middle 
particle representing the center of mass of the dot. This spectrum transforms to that of a "floppy" molecule 
for smaller value of Rw (i.e., for stronger confinements characterized by a larger value of ujo, and/or for 
weaker inter-electron repulsion), ultimately converging to the independent-particle picture associated with 
the circular central mean-field of the QD. Further evidence for the formation of the electron molecule was 
found through examination of the conditional probability distribution (CPD); that is, the anisotropic pair 
correlation P(r,ro), which expresses the probability of finding a particle at r given that the "observer" 
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(reference point) is riding on another particle at Tq EEl, 

N N 

P(r, ro) = (*| ^ ^ ^(r. - r)6{r, - ro)m/{^\^). (1) 

Here vjf (n, r2, ..., rjv) denotes the many-body wave function under consideration. It is instructive to note 
here certain similarities between the formation of a "two-electron molecule" in man-made quantum dots, 
and the collective (rovibrational) features observed in the electronic spectrum of doubly-excited heUum 
atoms ilEIIIl. 

For confined (finite) systems with a larger number of particles, one must resort to approximate computa- 
tional schemes. Of particular interest are methodologies that permit systematic evaluation of high-accuracy 
solutions to these many -body strongly-correlated systems, under field-free conditions, as well as under the 
influence of an applied magnetic field. We remark here, that the relatively large (spatial) size of QDs (re- 
sulting from materials' characteristics, e.g., a small electron effective mass and large dielectric constant), 
allows the full range of orbital magnetic effects to be covered for magnetic fields that are readily attained 
in the laboratory (less then 40 T). In contrast, for natural atoms and molecules, magnetic fields of sufficient 
strength (i.e., larger than 10^ T) to produce novel phenomena related to orbital magnetism (beyond the per- 
turbative regime), are known to occur only in astrophysical environments (e.g., on the surface of neutron 
stars). 

The 2D hamiltonian of the problem under consideration is given by 

■^-^ 2m* V c / ^-^ 2 ^-^ ^-^ KVij 

i—l i—1 i—1 j>i -^ 

describing N electrons (interacting via a Coulomb repulsion) confined by a parabolic potential of frequency 
uiQ and subjected to a perpendicular magnetic field B, whose vector potential is given in the symmetric 
gauge by A(r) = B x r/2 = {—By,Bx,0)/2. For sufficiently high magnetic field values (i.e., in the 
fractional quantum Hall effect, or FQHE, regime), the electrons are fully spin-polarized and the Zeeman 
term (not shown here) does not need to be considered. In the B ^ oo limit, the external confinement ujq 
can be neglected, and H can be restricted to operate in the lowest Landau level (LLL), reducing to the form 

OEiiniiMi 



^-i^-^^ + EE^- (3) 

where Wc = eB/{m*c) is the cyclotron frequency. 

For finite N, the solutions to the Schrodinger equations corresponding to the hamiltonians given by 
EqQ (with or without a magnetic field), or by Eq(|3} (in the i? — > cx) limit), must have a good angular 
momentum, L, and good spin quantum numbers (the latter is guaranteed in the high B, fully spin-polarized 
case). As described in detail elsewhere llllll2l[T^ll4l [T?i 16 1, these solutions can be well approximated 
by a two-step method consisting of symmetry breaking at the spin-and-space unrestricted Hartree-Fock 
level and subsequent symmetry restoration via post-Hartree-Fock projection techniques. We recall that 
the sS-UHF method relaxes both the double-occupancy requirement (namely, different spatial orbitals 
are employed for different spin directions), as well as relaxing the requirement that the electron (spatial) 
orbitals be constrained by the symmetry of the external confining potential. 

Results obtained for various approximation levels for a two-electron QD with B = and Rw = 
2.40 (that is, in the Wigner-molecule regime) are displayed in Fig. 2. In these calculations 1 16 1, the spin 
projection was performed following reference L17J . i.e., one constructs the wave function 

«'Sp(s)=7'.pi„(s)*uHF, (4) 
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Fig. 2 Various approximation levels for a field-free two-electron QD with Rw ~ 2.40. (a) Electron density of 
the RHF solution, exhibiting circular symmetry (due to the imposed symmetry restriction). The correlation energy 
Ec = 2.94 meV, is defined as the difference between the energy of this state and the exact solution [shown in frame (e)]. 
(bl) and (b2) The two occupied orbitals (modulus square) of the symmetry-broken "singlet" sS-UHF solution (bl), 
with the corresponding total electron density exhibiting non-circular shape (b2). The energy of the sS-UHF solution 
shows a gain of 44.3% of the correlation energy, (c) Electron density of the spin-projected (SP) singlet, showing broken 
spatial symmetry, but with an additional gain of correlation energy, (d) the spin-and-angular-momentum projected state 
(S&AMP) exhibiting restored circular symmetry with a 73.1% gain of the correlation energy. The choice of parameters 
is: dielectric constant n — S, parabolic confinement ujo = 5 meV, and effective mass m* — 0.067TTie. Distances are 



where ^uhf is the original symmetry-broken UHF determinant. In Eq. (|4}, the spin projection operator is 
given by 

s^ - s'{s' + i)n^ 



'Pspinis) = n 



s'^s 



[sis + 1) - s'{s' + l)]h 



2 ' 



(5) 



where the index s' runs over the quantum numbers of S"^, with S being the total spin. 
The angular momentum projector is given by 

p2lT 



(6) 



where L = /i + ^2 is the total angular momentum operator. As seen from Eq. (|6j, application of the projec- 
tion operator Vl to the spin-restored state ^sp(s) corresponds to a continuous configuration interaction 
(CCI) formalism. 

In the following we focus on the ground state of the system, i.e., L = 0. The energy of the projected 
state is given by 

£;pRj {L) = y "" h{-i)e''"'d-il I " n(7)e'^^d7, (7) 

with /i (7) = (*sp(s;0)|i/|*sp(s;7)) and 71(7) = (*sp(s; 0)|*sp(s; 7)), where *sp(s; 7) is the spin- 
restored wave function rotated by an azimuthal angle 7 and H is the many-body hamiltonian. We note that 
the UHF energies are simply given by -Euhf ~ h{0)/n{0). 

The electron densities corresponding to the initial RHF approximation [shown in Fig. 2(a)] and the 
final spin-and-angular-momentum projection (S&AMP) [shown in Fig. 2(d)], are circularly symmetric, 
while those corresponding to the two intermediate approximations, i.e., the sS-UHF and spin-projected 
(SP) solutions [Figs. 2(b2) and 2(c), respectively] break the circular symmetry. This behavior illustrates 
graphically the meaning of the term restoration of symmetry, and the interpretation that the sS-UHF broken- 
symmetry solution refers to the intrinsic (rotating) frame of reference of the electron molecule. In light of 
this discussion the final projected state is called a rotating Wigner molecule, or RWM. 



© 2003 WILEY- VCH Veikg GmbH & Co. KGaA, Wcinheim 



pss header will be provided by the publisher 



RWM 



EXD 



RWM 



EXD 



in 







■^^^ outer -^^Jf^ 



o 
II 




N=9 




Fig. 3 Conditional probability distributions (CPDs) at high B, evaluated for parabolic quantum dots through: (i) the 
two-step procedure of symmetry breaking and subsequent restoration, resulting in rotating Wigner molecules (RWM) 
(shown in the left columns for given A'^), and (ii) exact diagonalization (EXD, shown in the right columns for given 
A'^). The angular momentum values, and corresponding values of the fractional filling [see Eqjsj], are given on the left. 
The optimal polygonal structure for a given A*' is given by (ni, 712) with rzi + 712 = N. For TV = 7, 8, and 9, these 
arrangements are (1,6), (1,7), and (2,7), respectively. The reference point for the calculation of the CPD is denoted by 
a filled dot. Note in particular the two CPDs shown for A'^ = 9, illustrating that for a reference point located on the 
outer ring, the inner ring appears uniform, and vice versa for a reference point located on the inner ring (bottom row 
on the right). These results illustrate that the rings rotate independently of each other. 



1 .2 Electrons in circular quantum dots under high magnetic fields 

To illustrate the emergence of RWMs in parabolically confined QDs under high B, we show in Fig. 3 
results obtained |14| through the aforementioned two-step computational technique for iV = 7, 8, and 9 
electrons, and compare them with the results derived from exact diagonalization of the Hamiltonian [see 
Eq(|2|i]. Systematic investigations of QDs under high B revealed electronic states of crystalline character. 
These states are found for particular "magic" angular momentum values (L) that exhibit enhanced stability 
and are called cusp states. For a given value of B, one of these L's corresponds to the global minimum, 
i.e., the ground state, and varying B causes the ground state and its angular momentum to change. The 
cusp states have been long recognized 1 181 as the finite- A^ precursors of the fractional quantum Hall states 
in extended systems. In particular, the fractional fillings i' (defined in the thermodynamic limit) are related 
to the magic angular momenta of the finite-iV system as follows 1191 



(8) 



2i ■ 



In the literature of the fractional quantum Hall effect (FQHE), ever since the celebrated paper 1201 by 
Laughlin in 1983, the cusp states have been considered to be the antithesis of the Wigner crystal and to 
be described accurately by liquid-like wave functions, such as the Jastrow- Laughlin (JL) |20 21 J and 
composite-fermion (CF) 1221 1231 ones. This view, however, has been challenged recently 1111 1121 by 
the explicit derivation of trial wave functions for the cusp states that are associated with a rotating Wigner 
molecule. As discussed elsewhere 1 11 12 1, these parameter-free wave functions, which are by construction 
crystalline in character, have been shown to provide a simpler and improved description of the cusp states, 
in particular for high angular momenta (corresponding to low fractional fillings). 

The crystalline arrangements that were found 1 1 1 , :il2, 13 14| consist of concentric polygonal rings 
[see the conditional probabilities displayed in Fig. 3, with the reference (observation) point denoted by 
a filled dot]. These rings rotate independently of each other (see in particular the two cases shown for 
N — 9), with the electrons on each ring rotating coherently 1 14| . The rotations stabilize the RWM relative 
to the static one - namely, the projected (symmetry restored) states are lower in energy compared to the 
broken-symmetry ones (the unrestricted Hartree-Fock solutions). 
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1 .3 Symmetry breaking of trapped bosons 

In closing this section, we remark that the emergence of crystalline geometric arrangements, discussed 
above for electrons confined in 2D quantum dots, appears to be a general phenomenon that is predicted 
1241 to occur also for trapped bosonic atomic systems (neutral or charged) when the interatomic repulsion 
is tuned to exceed the characteristic energy of the harmonic trap. Indeed, application of the aforemen- 
tioned two-step method, allowed the evaluation 1 24l of solutions to the many -body Hamiltonians describing 
such bosonic systems, going beyond the mean-field (Gross-Pitaevskii) approach. These highly-correlated 
trapped 2D states exhibit localization of the bosons into polygonal-ringlike crystalline patterns, thus ex- 
tending earlier work that predicted localization of strongly repelling ID bosons (often referred to as the 
Tonk-Girardeau regime 125, 261 ). to a higher dimension. It is expected that these theoretical findings will 
be the subject of experimental explorations, using methodologies that have led recently to observations of 
localization transitions in ID boson systems 1271 1281 . 

2 Localization and entanglement in a two-electron elliptic quantum dot 

As discussed in the previous section, electron localization leading to formation of molecular-like structures 
[the aforementioned Wigner molecules] within a single circular two-dimensional (2D) quantum dot at zero 
magnetic field (B) has been theoretically predicted to occur l4ll5l 171 1291 1301 . as the strength of the e-e 
repulsive interaction relative to the zero-point energy increases. Of particular interest is a two-electron 
(2e) WM, in light of the current experimental effort |31 32 1 aiming at implementation of a spin-based 
l33l soUd-state quantum logic gate that employs two coupled one-electron QDs (double dot). 

Here, we present an exact diagonalization (EXD) and an approximate (generalized Heitler-London, 
GHL) microscopic treatment for two electrons in a single elliptic QD specified by the parameters of a 
recently investigated experimental device 1341 . While formation of Wigner molecules in circular QDs 
requires weak confinement, and thus large dots of lower densities (so that the interelectron repulsion dom- 
inates), we show that formation of such WMs is markedly enhanced in highly deformed (e.g., elliptic) dots 
due to their lower symmetry. The calculations provide a good description of the measured J{B) curve (the 
singlet-triplet splitting) when screening I35II36I due to the metal gates and leads is included (in addition to 
the weakening of the effective inter-electron repulsion due to the dielectric constant of the semiconductor, 
GaAs). In particular, our results reproduce the salient experimental findings pertaining to the vanishing of 
J{B) for a finite value of i? ~ 1.3 T [associated with a change in sign of J{B) indicating a singlet-triplet 
(ST) transition], as well as the flattening of the J{B) curve after the ST crossing. These properties, and in 
particular the latter one, are related directly to the formation of an electron molecular dimer and its effec- 
tive dissociation for large magnetic fields. The effective dissociation of the electron dimer is most naturally 
described through the GHL approximation, and it is fully supported by the more accurate, but physically 
less transparent, EXD. 

Of special interest for quantum computing is the degree of entanglement exhibited by the two-electron 
molecule in its singlet state |33J. Here, in relation to the microscopic calculations, we investigate two 
different measures of entanglement 1371 . The first, known as the concurrence (C) for two indistinguishable 
fermions [38"391, has been used in the analysis of the experiment in Ref. [34] (this measure is related to the 
operational cycle of a two-spin-qubit quantum logic gate 138 391 ). The second measure, referred to as the 
von Neumann entropy (S) for indistinguishable particles, has been developed in Ref. 1401 and used in Ref. 
1411. We show that the present wave-function-based methods, in conjunction with the knowledge of the dot 
shape and the J{B) curve, enable theoretical determination of the degree of entanglement, in particular for 
the elliptic QD of Ref. 1341 . The increase in the degree of entanglement (for both measures) with stronger 
magnetic fields correlates with the dissociation of the 2e molecule. This supports the experimental assertion 
1341 that cotunneling spectroscopy can probe properties of the electronic wave function of the QD, and not 
merely its low-energy spectrum. Our methodology can be straightforwardly applied to other cases of 
strongly-interacting devices, e.g., double dots with strong interdot-tunneling. 
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2.1 Microscopic treatment 
The Hamiltonian for two 2D interacting electrons is [see Eq. (|2j] 

n = H{Yi)+H{v2)+e'/{Kri2), (9) 

where the last term is the Coulomb repulsion, k, is the dielectric constant, and ri2 = |ri — r2|. H{y) 
is the single-particle Hamiltonian for an electron in an external perpendicular magnetic field B and an 
appropriate confinement potential. When position-dependent screening is included, the last term in Eq. Q 
is modified by a function of ri2 (see below). For an elliptic QD, the single-particle Hamiltonian is written 

as 

H{v) =T+ \m*{ulx^ + uy) + ^B . s, (10) 

where T = (p — eA/c)^/2m*, with A = 0.5{— By, Bx, 0) being the vector potential in the symmetric 
gauge, m* is the effective mass and p is the linear momentum of the electron. The second term is the 
external confining potential; the last term is the Zeeman interaction with g* being the effective g factor, fiB 
the Bohr magneton, and s the spin of an individual electron. 

The GHL method for solving the Hamiltoninian (|9jl consists of two steps. In the first step, we solve 
selfconsistently the ensuing unrestricted Hartree-Fock (UHF) equations allowing for lifting of the double- 
occupancy requirement (imposing this requirement gives the restricted HF method, RHF). For the Sz = 
solution, this step produces two single-electron orbitals UL,fl(r) that are localized left (L) and right (i?) 
of the center of the QD [unlike the RHF method that gives a single doubly-occupied elliptic (and sym- 
metric about the origin) orbital]. At this step, the many-body wave function is a single Slater determi- 
nant ^uhf(1 Tj2 i) = |ml(1 '\)uii(2 D) made out of the two occupied UHF spin-orbitals ul{1 t) = 
UL{ri)a{l) and u/j(2 |) = uji{r2)(3{2), where a{(3) denotes the up (down) ["f (J,)] spin. This UHF 
determinant is an eigenfunction of the projection Sz of the total spin S — si + S2, but not of S'^ (or the 
parity space-reflection operator). 

In the second step, we restore the broken parity and total-spin symmetries by applying to the UHF 
determinant the projection operator 1 15 42 1 V^'^^ = 1 T ^12, where the operator zui2 interchanges the 
spins of the two electrons [this is a special case of the operator given in Eq. (|5|i]; the upper (minus) sign 
corresponds to the singlet. The final result is a generalized Heitler-London (GHL) two-electron wave 
function ^qjjl(i"1) i"2) for the ground-state singlet (index s) and first-excited triplet (index t), which uses 
the UHF localized orbitals, 

*c'HL(i-i,r2) ex {uL{ri)uii{r2)±ULir2)uB.{ri))x"'\ (H) 

where x*'* = ('^(l)/3(2) T Q;(2)/3(l)) is the spin function for the 2e singlet and triplet states. The general 
formalism of the 2D UHF equations and of the subsequent restoration of broken spin symmetries can be 
found in Refs. lT5llT6ll29!'42l. 

The use of optimized UHF orbitals in the GHL is suitable for treating single elongated QDs. The GHL 
is equally applicable to double QDs with arbitrary interdot-tunneling coupling 1151 1421 . In contrast, the 
Heitler-London (HL) treatment 1 43 1 (known also as Valence bond), where non-optimized "atomic" orbitals 
of two isolated QDs are used, is appropriate only for the case of a double dot with small interdot-tunneling 
couphng i^ . 

The orbitals m^ fl;(r) are expanded in a real Cartesian harmonic-oscillator basis, i.e., 

K 

UL,fl(r) = 5]cf«^,(r), (12) 

where the index j = (m, n) and (pj{r) — Xm{x)Yn{y), with Xm{Yn) being the eigenfunctions of the 
one-dimensional oscillator in the x(y) direction with frequency LOxipJy). The parity operator V yields 
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Fig. 4 The singlet-triplet splitting J = E" — E* as z function of the magnetic field B for an elliptic QD with 
huix = 1.2 meV and huy — 3.3 meV (these values correspond to the device of Ref. |34|). Solid line: GHL (broken- 
symmetry UHF -I- restoration of symmetries) results with a coordinate-independent screening (k = 22). Dashed line: 
EXD results with k = 12.9 (GaAs), but including screening with a coordinate dependence according to Ref. 1 36 1 and 
d = 18.0 nm (see text). The rest of the material parameters used are: m,*(GaAs)= 0.067me, and g* — (see text). 
The experimental measurements [34| are denoted by open squares. Our sign convention for J is opposite to that in 
Ref. Illl. 

VXmix) = {—l)"^Xm{x), and similarly for Yn{y). The expansion coefficients C ' are real for B = 
and complex for finite B. In the calculations we use K = 79, yielding convergent results. 

In the EXD method, the many-body wave function is written as a linear superposition over the basis of 
non-interacting two-electron determinants, i.e.. 



*EXD(l"l,r2) 



2K 

E 



n^^f\^P{l■,^M2■,J)), 



(13) 



where ^p{l;i) = ifi{l t) if 1 < « < if and V(l;«) = ^i-k{1 i) if K + 1 < i < 2K [and similarly 
for ?/'(2, j)]. The total energies i?|xD '^^'^ ^^^ coefficients fi^-'- are obtained through a "brute force" diago- 
nalization of the matrix eigenvalue equation corresponding to the Hamiltonian in Eq. (|9}. The EXD wave 
function does not immediately reveal any particular form, although, our calculations below show that it 
can be approximated by a GHL wave function in the case of the elliptic dot under consideration. 

To model the experimental elliptic QD device, we take, following Ref. [34], hujx = 1.2 meV and hujy = 
3.3 meV. The effective mass of the electron is taken as m* = Q.OGTnie (GaAs). Since the experiment did 
not resolve the lifting of the triplet degeneracy caused by the Zeeman term, we take g* = 0. Using the two 
step method, we calculate the GHL singlet-triplet splitting Jghl(-B) = Eqj^i^{B) — EQj^f^{B) as a function 
of the magnetic field in the range < B < 2.5 T. Screening of the e-e interaction due to the metal gates 
and leads must be considered in order to reproduce the experimental J{B) curve 1441 . This screening can 
be modeled, to first approximation, by a position-independent adjustment of the dielectric constant k 1451 . 
Indeed, with k — 22.0 (instead of the GaAs dielectric constant, i.e., k — 12.9), good agreement with the 
experimental data is obtained [see Fig. 4]. In particular, we note the singlet-triplet crossing for B Ki 1.3 T, 
and the flattening of the J{B) curve beyond this crossing. 

We have also explored, particularly in the context of the EXD treatment, a position-dependent screening 
using the functional form, {e^ / Kri2)[l — (1 + 4d^/r^2)~^^^]' proposed in Ref. |36|, with d as a fitting 
parameter. The Jexd{B) result for d — 18.0 nm is depicted in Fig. 4 (dashed line), and it is in very good 
agreement with the experimental measurement. 

The singlet state electron densities from the GHL and the EXD treatments at i? = and B — 2.5 T 
are displayed in Fig. 5. These densities illustrate the dissociation of the electron dimer with increasing 
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GHL B=0 




GHL B=2.5T 





Fig. 5 Total electron densities (EDs) associated with the singlet state of the elliptic dot at B = and _B = 2.5 T. (a) 
The GHL densities, (b) The EXD densities. The rest of the parameters and the screening of the Coulomb interaction 
are as in Fig. 4. Lengths in nm and densities in 10^'* nm^^. 



magnetic field. The asymptotic convergence (beyond the ST point) of the energies of the singlet and triplet 
states, i.e., [J{B) — > as i? ^ cxj] is a reflection of the dissociation of the 2e molecule, since the ground- 
state energy of two fully spatially separated electrons (zero overlap) does not depend on the total spin 

ESI. 



2.2 Measures of entanglement 1371 

To calculate the concurrence C 138* '391, one needs a decomposition of the GHL wave function into a 
linear superposition of orthogonal Slater determinants. Thus one needs to expand the nonorthogonal 
u^'^{t) orbitals as a superposition of two other orthogonal ones. To this effect, we write M^'^(r) ex 
'I'"'" (r) ± ^$^ (r), where <!>+ (r) and $^ (r) are the parity symmetric and antisymmetric (along the x-axis) 
components in their expansion given by Eq. \\2l . Subsequently, with the use of Eq. Jilt , the GHL singlet 
can be rearranged as follows: 



4-, 



GHL 



(x|$+(lt)$+(2i)) -,;!$-(! t)<i>-(2i)), 



(14) 



where the so-called interaction parameter II39I ,77 = ^^ , is the coefficient in front of the second determinant. 
Knowing r] allows a direct evaluation of the concurrence of the singlet state, since C" = 277/(1 + rf^) 
|39|. Note that ^^(r) and $^(r) are properly normalized. It is straightforward to show that 77 = (1 — 
\Slr\)I{'^ + \Slr\), where Slr (with 15^^] < 1) is the overlap of the original u^'^{r) orbitals. 
For the GHL triplet, one obtains an expression independent of the interaction parameter 77, i.e.. 



*, 



GHL 



cx|<l>+(lt)<i>-(2i)) + |$+(U)<l>-(2t)), 



(15) 



which is a maximally (C* — 1) entangled state. Note that underlying the analysis of the experiments in Ref. 
1341 is a conjecture that wave functions of the form given in Eqs. (I14> and (I15> describe the two electrons 
in the elliptic QD. 

For the GHL singlet, using the overlaps of the left and right orbitals, we find that starting with rj — 0.46 
(C = 0.76) at i? = 0, the interaction parameter (singlet-state concurrence) increases monotonically to 
7] = 0.65 (C — 0.92) at _B = 2.5 T. At the intermediate value corresponding to the ST transition {B ^ 1.3 
T), we find -q — 0.54 (C* ~ 0.83) [47]. Our B — theoretical results for 77 and C* are in remarkable 
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Fig. 6 Von Neumann entropy for the singlet state of the elliptic dot as a function of the magnetic field B. Solid line: 
GHL. Dashed line: EXD. The rest of the parameters and the screening of the Coulomb interaction are as in Fig. 4. On 
the left, we show histograms for the \zu^ coefficients [see Eq. <16H of the singlet state at B = 1.3 T, illustrating the 
dominance of two configurations. Note the small third coefficient \z-i\^ = 0.023 in the EXD case. 



agreement with the experimental estimates L34J of 77 = 0.5 ± 0.1 and C = 0.8, which were based solely 
on conductance measurements below the ST transition (i.e., near B = 0). 

To compute the von Neumann entropy, one needs to bring both the EXD and the GHL wave functions 
into a diagonal form (the socalled "canonical form" I40. .48J ). i.e.. 



M 



*EXD(ri, r2) = Y. ^fe'*l'^(l; 2fc - l)<i>(2; 2fc)), 



(16) 



k=\ 



with the $(i)'s being appropriate spin orbitals resulting from a unitary transformation of the basis spin 
orbitals ip{jys [see Eq. ( I13H : only terms with Zk ^ contribute. The upper bound M can be smaller 
(but not larger) than K (the dimension of the single-particle basis); M is referred to as the Slater rank. 
One obtains the coefficients of the canonical expansion from the fact that the |zfcp are eigenvalues of 
the hermitian matrix ft'Tl [J7, see Eq. (I13> . is antisymmetric]. The von Neumann entropy is given by 



S 






^M 



Zk\ \og2i\zk\ ) with the normalization ^^_j^ \zk\ — 1- Note that the GHL wave functions 
in Eqs. ( I14> and ( I15> are already in canonical form, which shows that they always have a Slater rank of 
M = 2. One finds Sp,„, = log2(l + ff) - rj'^ log2(?7^)/(l + r/^), and 5^jjl = 1 for ^H B. For large B, 



GHL 

the overlap between the two electrons of the dissociated dimer vanishes, and thus 77 ^ 1 and 5qjjl ^ 1- 

Since the EXD singlet has obviously a Slater rank M > 2, the definition of concurrence is not applicable 
to it. The von Neumann entropy for the EXD singlet (5|xd) is displayed in Fig. 6, along with that (Sq^^^) of 
the GHL singlet. iS|xd ^"^d iSqjjl are rather close to each other for the entire B range, and it is remarkable 
that both remain close to unity for large B, although the maximum allowed mathematical value is logj (K) 
[as aforementioned we use A' = 79, i.e., log2(79) — 6.3]; this maximal value applies for both the EXD and 
GHL approaches. The saturation of the entropy for large B to a value close to unity reflects the dominant 
(and roughly equal at large B) weight of two configurations in the canonical expansion [see Eq. ( I16H of 
the EXD wave function, which are related to the two terms (M ~ 2) in the canonical expansion of the 
GHL singlet [Eq. ( I14H . This is illustrated by the histograms of the \zf.\'^ coefficients for B = 1.3 T in 
Fig. 6 (left column). These observations support the GHL approximation, which is computationally less 
demanding than the exact diagonalization, and can be used easily for larger N. 

3 Summary 

We discussed symmetry breaking in two-dimensional quantum dots resulting from strong interelectron 
repulsion relative to the zero-point kinetic energy associated with the confining potential. Such symmetry 
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breaking leads to the emergence of crystalline arrangements of electrons in the dot. The so-called Wigner 
molecules form already at field-free conditions. The appearance of rotating Wigner molecules in circular 
dots under high magnetic field, and their relation to magic angular momenta and quantum-Hall-effect 
fractional fillings was also discussed. 

Furthermore, we have shown, through exact and approximate microscopic treatments, formation of an 
electron molecular dimer in an elliptic QD (Fig. 5) for screened interelectron repulsion. The formation 
and effective dissociation (in high magnetic fields) of the electron dimer are reflected in the behavior of 
the computed singlet-triplet splitting, J{B), that agrees well (see Fig. 4) with the measurements 1341 . 
Furthermore, we showed that, from a knowledge of the dot shape and of J{B), theoretical analysis along 
the lines introduced here allows probing of the correlated ground-state wave function and determination 
of its degree of entanglement. This presents an alternative to the experimental study where determination 
of the concurrence utilized conductance data L34J . We have employed two measures of entanglement for 
indistinguishable fermions (the concurrence and the von Neumann entropy) and have shown that their 
behavior correlates with the effective dissociation of the electron dimer. Such information is of interest to 
the implementation of spin-based solid-state quantum logic gates. 

This research is supported by the US D.O.E. (Grant No. FG05-86ER45234), and NSF (Grant No. 
DMR-0205328). We thank M. Pustilnik for comments on the manuscript. 
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